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Abstract 

We study single-electron quantum dots on helium surface created by electrodes submerged into the helium. The intradot 
potential is electrostatically controlled. We find the electron energy spectrum and identify relaxation mechanisms. Strong 
in-plane confinement significantly slows down electron relaxation. Energy relaxation is due primarily to coupling to phonons 
in helium. Dephasing is determined by thermally excited ripplons and by noise from underlying electrodes. The decay rate 
can be further suppressed by a magnetic field normal to the helium surface. Slow relaxation in combination with control 
over the energy spectrum make localized electrons appealing as potential qubits of a quantum computer. 
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The system of electrons on the surface of supcrfiuid 
4 He is attractive from the point of view of making a 
scalable quantum computer (QC) [1]. The electrons 
have an extremely long relaxation time and display the 
highest mobility known in a condensed-matter system 
[2] . A QC can be made by submerging a system of in- 
dividually addressed micro-electrodes beneath the he- 
lium surface. The typical interelectron distance is com- 
paratively large, ~ 1 /im, which simplifies fabrication 
of an electrode array [3]. The electrode potential, the 
high barrier that prevents electrons from penetrating 
into the helium, and the helium image potential to- 
gether create a single-electron quantum dot above each 
electrode. The parameters of the dot can be controlled 
by the electrode potential. 

Here we study the energy spectrum and dissipation 
processes for electrons in quantum dots on helium sur- 
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face. We discuss mechanisms of coupling to helium ex- 
citations, phonons and ripplons, and the dependence 
of the electron relaxation rate on the quantum dot pa- 
rameters. We investigate the role of a magnetic field 
normal to the surface and of the electron-electron in- 
teraction. Decay and decoherence of the electron states 
result also from classical and quantum electrode noise. 
We relate the corresponding relaxation rates to the 
power spectrum of the fluctuating electric field on the 
electron and analyze their dependence on the parame- 
ters of the electrodes and external leads. 

The geometry of a quantum dot can be understood 
from Fig. 1. The potential is a sum of the out-of-plane 
and in-plane parts. The out-of-plane potential is simi- 
lar to that in the absence of the electrode. It leads to 
quantization of motion normal to the surface. In the ab- 
sence of the field E± the energy levels are E„ = —R/n 2 
(n — 1,2,...), where the effective Rydberg energy 7? w 
8 K [4]. The in-plane potential Un(r) « muj 2 r 2 /2 [r = 
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(x, y)] is parabolic, for typical parameter values. 

In the absence of a magnetic field, electron states 
in a dot \n,v,m v ) are characterized by the quantum 
number n of motion normal to the surface, the principal 
quantum number v of vibrations about an equilibrium 
in-plane position, and the number m„ that enumerates 
degenerate vibrational states. 

A confined electron can serve as a qubit [1]. The 
working states of the qubit are 1 1, 0, 0} and |2, 0, 0) . The 
energy difference between these states Ei — E\ can be 
Stark-shifted by the electric field from the electrode 
£±. The shift of 1 GHz occurs if £± is changed by 
~ 1 V/cm. The field £± also determines the in- plane 
vibrational frequency w,, . A simple estimate of £± can 
be made by assuming that the electrodes are spheres 
of radius r e i. In this case u).. = {e£±/mh) 1/2 . Typically 
W|| /2n ~ 20 GHz, whereas the transition frequency 
f2tr = {E2 — Ei)/h is 6-10 times larger. 




Fig. 1. Geometry of a single-electron quantum dot on he- 
lium surface. A micro-electrode is submerged by the depth 
h ~ 0.S fim beneath the helium surface. The electron is 
driven by a field E± normal to the surface. This field is 
a sum of the electrode field £j_ and the field of the paral- 
lel-plate capacitor (only the lower plate of the capacitor is 
shown). The confining in-plane potential U» (r) [r = (x, y)\ 
is determined by the electrode potential and geometry. 

Because of the discreteness of the electron energy 
spectrum in a dot, the mechanisms of electron decay 
and dephasing are qualitatively different from those 
studied for a 2D electron system on helium in the ab- 
sence of in-plane confinement [5] . In particular, the ma- 
jor known scattering mechanism, quasi-elastic scatter- 
ing by ripplons, does not work. Decay of the excited 
state |2,0,0) is most likely to occur via a ripplon- or 
phonon-induced transition to the closest in energy ex- 
cited vibrational state of the electron |1, u c ,m v ), with 
v c — mt[(E2 — Ei)/huA. The energy transfer in the 



transition is SE — E2 — E\ — vJvjJ,, ~ hu,, . It largely 
exceeds the energy of ripplons with wave numbers q < 
I/ay , where a,, = (h/mu)» ) 1//2 is the in-plane electron 
localization length in the dot. This makes one-ripplon 
decay exponentially improbable and strongly reduces 
the decay rate compared to the case of electrons that 
are free to move along the helium surface. 

An electron in a dot can decay by emitting two 
ripplons that will propagate in opposite directions 
with nearly same wave numbers q\,2, which are deter- 
mined by the condition that the ripplon frequency is 
Wr(gi,2) ~ SE/2h [|qi + q 2 | < l/a {[ < <?i, 2 ]. Alterna- 
tively, and even with higher probability, decay may 
occur through an electron transition accompanied by 
emission of a phonon. The appropriate phonons prop- 
agate nearly normal to the helium surface: their wave 
vectors make an angle with the normal to the surface 
~ (mVg/hu},, ) 1//2 <C 1, where v s is the sound velocity 
in helium. We propose two mechanisms of coupling to 
phonons [6] . One is phonon-induced deformation of the 
helium surface. The second is phonon-induced modu- 
lation of the dielectric constant of helium and thus of 
the electrostatic energy of the electron above helium. 

Relatively simple expressions for the electron decay 
rate for different coupling mechanisms are obtained 
by approximating the electron potential at the helium 
surface by a sharp infinitely high potential wall. Such 
an approximation applies only if the wavelength of he- 
lium vibrations is much bigger than the width Ad of 
the diffuse surface layer, which is of the order of a few 
angstroms [7]. When we estimate the decay rate nu- 
merically, we make an assumption that helium vibra- 
tions are essentially decoupled from the electrons when 
their wave numbers exceed AT 1 . Then, for typical in- 
plane electron frequencies w,, /2-7T ~ 20 GHz, both the 
two-ripplon and phonon decay rates are determined 
primarily by electron transitions with the smallest en- 
ergy transfer 8E ~ hu),, . Transitions over several elec- 
tron vibrational levels, with energy transfer ntiw,, with 
n>l, can be disregarded. This gives the decay rate 
< 10 4 s _1 , which is presumably an overestimate. 

Electron dephasing due to coupling to excitations 
in helium comes primarily from quasi-elastic scatter- 
ing of ripplons off the electron. Scattering amplitudes 
are different in different electron states, and therefore 
scattering leads to diffusion of the phase difference be- 
tween the wave functions |2, 0, 0) and |1, 0, 0). The typ- 
ical wave numbers of the ripplons that contribute 
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to dephasing are determined by the condition that the 
ripplon frequency u) r {q^) < ksT/h. They are much 
smaller than 1/A d , and the coupling is well described 
by the approximation of a sharp helium surface. For dif- 
ferent mechanisms of coupling to ripplons the dephas- 
ing rate r'*' displays same temperature dependence 
T^' oc T 3 , which is much slower than the standard T 7 
law for defects in solids. Numerically, F^' < 10 2 s _1 
for T — 10 mK. 

The electron decay rate may be further reduced by 
applying a magnetic field normal to the helium surface. 
As we will show, this is equivalent to further increas- 
ing u>n and should require helium vibrations with even 
shorter wavelengths in order to meet the condition of 
energy conservation in decay. The theory will also de- 
scribe the case of many quantum dots. In this case, the 
energy spectrum of in-plane electron excitations con- 
sists of plasmon bands and is continuous. Still, as we 
show, the relaxation rate remains strongly suppressed 
compared to the case of unconfined electrons. 

The in-plane electron coordinate r n in an nth dot 
can be expanded in the creation and annihilation oper- 
ators of the in-plane vibrational modes of the electrons 

r n =R„ + ^^[Ag ) a fcj +H.c]. (1) 

Here, R n is the equilibrium in-plane position, and j — 
1, 2 enumerates vibrational modes. The quantum num- 
ber k can be set equal to zero in the case of one dot, 
whereas for a periodic array of dots it becomes a plas- 
mon wave vector k. 

For one dot (n = 1) in a strong magnetic field B 
antiparallel to z 

woi ~u) c = \eB\/mc, lu 02 ~ &o = w 2 /w c , (2) 
Ai « iA 2 » 2~ 1/2 /(x + iy) [A, = A<V], 

where I = (S/mw c )'' 2 ] is the magnetic length, and we 
assumed that Wy is small compared to the cyclotron 
frequency u> c . For a multi-dot system, the vibrational 
frequencies form two bands. The bandwidths are 
< LUp/uJc, where u) p = (2-7re 2 /md 3 ) 1 ' ' 2 is the typical 
plasma frequency {u) p -C cu c ) and d is the interelectron 
distance. They are further reduced if u) p <C ui^. The 
minimal frequency of the upper band is « ui c , whereas 
that of the lower band is ~ Qo- 

The Hamiltonian that describes 1 2) — > 1} transitions 
induced by excitations in helium has the form 



#} d) = E I 2 >«"UIE ^ <qr "+H.c. (3) 

Here, |1)„ and |2)„ are the states of an nth electron nor- 
mal to the surface, and Vq is the operator that depends 
on the coordinates of helium vibrations, i.e., phonons 
and ripplons. The wavelengths of the vibrations in- 
volved in electron scattering are much smaller than the 
interelectron distance. Therefore each electron has its 
"own" thermal bath of helium excitations. In the Born 
approximation, the decay rate of the state |2}„ for an 
nth electron is 

ri d) = ?T 2 Re / dte lSltrt J2 <5„(q,t) 

X<V q (t)VL q (0)>, <S Il (q,t) = (e lqr " (i) e- lqr "W>. (4) 

Here, the averaging (•) is performed assuming that the 
electron and helium vibrations are uncoupled. We have 
also disregarded the difference between the transition 
frequencies Q tT — (E-2 — E\)/h for different electrons. 

From Eq. (1), for ksT <§C hQo the structure factor 
<S n (q, i) (4) has the form 

<S„(q,t) =exp[-qV„(i)/2], (5) 
W n (t) = J2 k ,\A^f[l-eM-i^t)}. 

The effective Debye- Waller factor W„(t) is indepen- 
dent of the electron number n if the coefficients A^ 
for different electrons differ only by a phase factor, as 
in the case of a periodic set of dots. 

In what follows we will assume that the level spac- 
ing E-2 — Ei is of the same order as the distance be- 
tween the Landau levels fkj c and that the energy deficit 
8E = E2 — Ei — v c hu] c largely exceeds hflo [here, i/ c — 
int(i?2 — Ei)/hio c ]. Then many vibrations of the lower 
vibrational branch j = 2 are excited in a transition. 
This means that the structure factor <S n (q, t) can be 
evaluated assuming that f2oi <C 1. In doing this we will 
take into account the relations J2 k lA^'l 2 w I 2 and 

£ lA^I 2 ^ » l 2 fl n , n n = [d 2 UJdr 2 n ]/2muj c , (6) 

k 

where U.. is the overall in-plane potential of an nth elec- 
tron that includes the Coulomb energy of the electron- 
electron interaction, and the Laplacian is calculated at 
the equilibrium position. It is seen from Eq. (2) that 
Eq. (6) applies in the case of one dot. A proof for a 
multi-dot system will be discussed elsewhere. 
From Eq. (6) we obtain 
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<S„(q,*) ^exph(g 2 / 2 /2)(l + ^nt)] 

^^(W^' (7) 

This equation shows that, in an electron |2) — > |1) 
transition, the energy transferred to the lower-branch 
vibrational modes is ~ q 2 l 2 Q. n h. The typical values of 
the factor q 2 l 2 are < 1, otherwise the transition prob- 
ability becomes exponentially small. This means that, 
even though the energy spectrum of a set of dots is 
band-like, the typical energy that has to be transferred 
to helium excitations is ~ 8E ~ f\io c . For u>^ <C uj c this 
is a much bigger energy than in the absence of a mag- 
netic field. Therefore we expect that a magnetic field 
can significantly reduce the decay rate, because decay 
requires helium excitations with wavelengths smaller 
than the width of the surface diffuse layer. 

The occurrence of a mode with frequency Qo -C w,, 
in a magnetic field leads to an increase of the ripplon- 
induced dephasing rate. This can be seen in the higher 
order in the electron-ripplon coupling where account 
is taken of virtual ripplon-induced transitions between 
the electron vibrational levels, cf. [1,6]. The amplitude 
of such transitions increases with the decreasing fio, 
i.e., with the increasing magnetic field. The effect im- 
poses a limitation on the field magnitude. One con- 
cludes that there exists an optimal range of magnetic 
fields where both decay and dephasing rates are small 
and of the same order of magnitude. 

An important source of dephasing of electron states 
in a single-electron dot is noise from the underlying 
electrode, see Fig. 1. We will discuss it for one dot, as 
the role of inter-dot interaction is minor in this case. 
Coupling to the electrode noise is dipolar, with the 
Hamiltonian 

-ftdip = eS£±z, (8) 

where S£± is the normal to the surface component of 
the field from quantum or classical charge density fluc- 
tuations on the electrode. 

The rates r'f and of the electron dephasing 
and decay can be expressed [6] in terms of the field 
correlation function 

roc 

Q{uj)= dte iwt (S£±{t)5£ ± (0)} (9) 
Jo 

as 

T ( f =e 2 (2 2 2-2ii) 2 ReQ(0)/ft 2 , 

=e 2 \z 12 \ 2 RoQ(Q tr )/h 2 . (10) 



Here, z i3 = {i\z\j), with i,j = 1,2. 

A major contribution to the dephasing rate comes 
from Johnson's noise in the external leads. A simple 
estimate can be made by assuming that the confining 
electrode is a sphere of radius r e i . It gives 

r<f> = 2fc B T cxt ft cxt e 2 (£ 22 - Zll ) 2 r 2 cl /h 2 h\ (11) 

where T cx t and 7?. C xt are the lead temperature and re- 
sistance. For = 25 Q, T ext = 1 K, r e i = 0.1 /tm, 
and h = 0.5 ^m we obtain T^f « 1 x 10 4 s _1 . Eq. (11) 
suggests how to reduce the rate F^'. The decay rate 
is much less than the phonon-induced decay rate. 

Eq. (10) allows one also to estimate the effect of 
noise from reorienting defects in the electrode. For 
electrodes submerged into helium this noise should be 
weaker than in semiconductor-based systems proposed 
for quantum computing, in particular because it scales 
with the distance to the electrode as h~ 4 , cf. Eq. (11). 

In this paper we have calculated the relaxation rate 
for electrons in quantum dots on the helium surface. 
We proposed new relaxation mechanisms and found 
the dependence of the relaxation rate on the potential 
of the confining electrodes and the magnetic field. It 
follows from the results that the ratio of the relaxation 
rate to the clock frequency of the QC based on electrons 
on helium, which is determined by the electron-electron 
interaction, can be as small as 10~ 4 — 10 -5 , for typical 
interdot spacing < 1 /urn. 
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